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d • Abstract 

O^. Double-well optical lattices are considered, each cite of which is formed by a double- 

well potential. The lattice is assumed to be in an insulating state and order and disorder 
are defined with respect to the displacement of atoms inside the double-well potential. 
It is shown that in such lattices, in addition to purely ordered and disordered states, 
there can exist an intermediate mixed state, where, inside a generally ordered lattice. 
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1 Introduction 

Cold atoms, loaded into optical lattices, form the systems that are unique with regard to 
the controllability of their properties (see the review articles [1-5] and references therein). 
One can vary the sorts of atoms as well as the lattice types and lattice spacing. For the 
given kind of atoms, one can vary their interactions by means of the Feshbach resonance 
techniques [5-7]. It is also possible to change the lattice properties by imposing external 
fields or shaking the lattice. This is why, one is able to create a variety of different system 
states in the lattice [1-5]. 

A special type of lattices are the double-well lattices, whose experimental realization has 
become possible in recent years [8-12]. In such lattices, each lattice site is formed by a double- 
well potential, which makes it feasible to create novel states that would be impossible in the 
standard single-well lattices. The phase diagram of bosons in double-well lattices has been 
studied [13,14], with the emphasis on the boundary of the insulator- sup erfiuid transition. 
Double-well lattices in an insulator state have been analyzed [15-17], with the emphasis on 
collective excitations in these lattices [15], the order- disorder transition [15], dynamics of 
atomic loading in the lattice [15,16], and the feasibility of regulating atomic imbalance by 
varying the double-well potential shape [17]. The possibility of transferring atoms between 
the wells of a double-well potential by means of the stimulated Raman adiabatic passage 
scheme has been suggested [18]. 

In the present paper, the possibility of creating another state in a double-well optical 
lattice is advanced. We consider an insulating double-well lattice, with a single atom per 
cite, which is a setup that would be typical for the Mott insulator in a single-well lattice. 
Because of the insulating state and the unity filling factor, the atom statistics are not crucial; 
the atoms can be bosons as well as fermions. In such a system, there can exist an ordered 
state, when all atoms are unidirectionally shifted into one of the wells of each of the double- 
well potentials [15]. Such a state is schematically shown in Fig. 1. It may also happen 
a disordered state, when atoms are randomly located in one of the wells of the double- 
well potential [15]. We analyze under what conditions there could arise an intermediate 
mixed state, when inside a generally ordered lattice there could appear mesoscopic regions of 
disordered states, as is illustrated in Fig. 2. The term mesoscopic leieis to the characteristic 
size Imes of the disordered regions, as compared to the lattice spacing a and to the total 
system size L. Of course, the disordered regions can be of different sizes and shapes, but it 
is always possible to define their characteristic, or average, size Imes- This size is mesoscopic 
in the sense of being intermediate between the lattice spacing and system sizes: 

CL <C l-rnes ^ ^ ■ 

We show that there exist such system parameters, when the intermediate mixed state of 
a lattice with mesoscopic disorder can arise. Since optical lattices allow for extremely wide 
possibilities in varying the system parameters, the suggested mixed state can be realized 
experimentally. And, as far as the properties of the mixed state can be regulated, it can be 
used for a variety of applications. 



2 Random Region Configurations 

In order to understand whether and when the mixed state could appear, it is necessary to 
take into account its possible existence, so that the probability of this state could be defined 
self-consistently. A system, formed by subregions corresponding to different thermodynamic 
states, having different orders, can be described by using the local-equilibrium Gibbs ensem- 
bles [19,20]. The regions of competing thermodynamics phases can arise in a self-organized 
way in arbitrary spatial locations, that is, their locations are randomly distributed in space. 
Therefore, it is necessary to average over these random configurations. Such a procedure 
has been suggested in Refs. [21-27] and reviewed in Ref. [28]. All mathematical details of 
this approach has been thoroughly described in these references. To better understand the 
following, in this section, we briefiy delineate the main points of the approach as applied to 
a lattice system. 

Let us consider the lattice 

Z = {a, : I = 1,2,..., Nl] (1) 

composed of A'"^ lattice cites characterized by the vectors a^. Suppose that there can exist 
two different thermodynamic states corresponding to an ordered and disordered phases. 
Meanwhile, it is not essential what kind of order is assumed. But for what follows, we keep 
in mind that each lattice cite is formed by a double-well potential and that the order is 
understood in the sense of the existence of an average atomic imbalance, as is explained in 
the Introduction. 

When a part of the lattice, say Zi, is filled by an ordered state and a part Z2, by a 
disordered one, the whole lattice is represented as the union 

Z = Zi|JZ2 (Zip|Z2 = o) , (2) 

in which 

Nl = N^ + N2 (N, = mes Z,) . 

This is what is called the orthogonal covering. Each manifold covering can be characterized 
[29] by the manifold characteristic functions, or manifold indicators 

^^'^ - \ , a, ^ Z, , ^^> 

describing the lattice separation into sublattices. The collection of the given manifold indi- 
cators composes the covering configuration 

e = {6.: 2 = 1, 2,..., AT,; z/ = l,2}. (4) 

Each thermodynamic state is characterized by the corresponding space of typical mi- 
crostates Tiu that is a weighted Hilbert space, the complete metric space with a weighted 
norm induced by the inner product [28]. The latter implies that, if {(/?„} is an orthonormal 
basis in l-Ly and / G l-iy, so that the inner product {(p, f) be given, then the weighted norm 
of / over Til, is 



'^Pnuli^nJW 



with the set {pni/} composing a probabiUty measure: 

<p„^ < 1 , "^Pnu = 1 • 

n 

The statistical average, with a statistical operator p^,, of an operator A on Tiu is defined as 

{A)^ = TlnJuA = ^Pnu ((fn^PuAipn) ■ 
n 

The ordered and disordered states are distinguished as follows. There should exist an 
order operator fj such that the averages 

Vv = (jl)v (5) 

define the order parameters for the corresponding ordered {rji ^ 0) and disordered {ri2 = 0) 
states. 

For each fixed configuration C,, the Hamiltonian energy operator, containing the single- 
atom terms Hi and the binary interaction terms Hij, possesses the representations 

defined on the related spaces "Hj,. The system as a whole is described by the Hamiltonian 

m) = m)®m) (7) 

defined on the mixture space 

M=Hi6dH2. (8) 



The regions of disorder can arise in arbitrary locations and have different shapes. There- 
fore, it is necessary to invoke an averaging procedure over the random configuration set (4). 
Denoting the stochastic averaging by the double angle brackets << ■ ■ ■ >>, one has the free 
energy 

F = -Tln((Tre~^^(«))) , (9) 

where the trace is over space (8) and f3 = 1/T is inverse temperature. This corresponds to 
annealed disorder, since the disorder regions are not fixed in space by an external force, but 
appear in the system in a self-organized way. The stochastic averaging of an operator A{C,) 
over the random configurations C, can be rigorously defined [28] as a functional integration 
over the manifold indicators, 

{{Mo))= I mvi- (10) 



Physically, the procedure of averaging over the random configurations is based on the ex- 
istence of different spatial scales related to atoms (microscopic scale) and to the regions of 
disorder (mesoscopic scale). Different spatial scales also are often connected with different 
temporal scales of slow and fast motion [30-33]. Technically, the averaging over configura- 
tions reminds the renormalization procedure [34,35]. All exact mathematical details of the 
configuration averaging is given in the review article [28]. 



Introducing the effective renormalized Hamiltonian 

H=-T\n{{e'^"^^^)) (11) 

makes it possible to rewrite the free energy (9) as 

F = -TlnTre"'^^ . (12) 

Thus, the basic point is to find the effective Hamihonian (11). For the Hamiltonian 
representation (6), there has been proved [27,28] the following statement. 

Theorem. The effective Hamiltonian (11), where H{C,) is given by the sum (7) of the 
Hamiltonian representations (6), takes the form 

^ = ^10^2, (13) 

in which 

H, = w,J2Hi + wlJ2 H^j , (14) 

with the weights w^ being the minimizers of the free energy (12), under the conditions 

wi + W2 = l , < w^ < 1 . (15) 

The Hamiltonian form (13) is similar to the channel Hamiltonian representation in scat- 
tering theory [36,37]. 

3 Double- Well Insulating Lattice 

The detailed derivation of the Hamiltonian for a double- well insulating lattice has been done 
in Ref. [15]. It is convenient to resort to pseudospin representation. To remind the reader 
the physical meaning of the pseudospin operators to be used in what follows, let us recall 
how they are constructed. 

We start with the operators of destruction and creation, c„j and c^^-, associated with the 
n-th energy level En of an atom in the j-th lattice cite. Considering an insulating lattice 
implies that atomic transitions between different lattice cites are suppressed, so that 

A lattice with the unity filling factor is assumed, when the number of atoms A^ coincides 
with the number of lattice cites N^. Then, limiting the consideration by two lowest energy 
levels of an atom in a double-well potential, one has the unipolarity condition 

^Ij^lj + '^2j'^2j = 1 5 CnjCnj = . (17) 



The pseudospin operators are introduced through the transformations 
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c[jC2j = S^ - iS] , c\fi, = S] + iS] . (18) 

The pseudospin operators 3°' satisfy the standard spin commutation relations for both the 
cases, when the operators c„j are either bosonic or fermionic. From the atomic operators, 
associated with the energy levels, one can pass to the operators, associated with the left and 
right positions of an atom in a double-well potential, by means of the transformations 

cii = -T^^^JL + Cjr) , C2j = -ds^jL - Cjr) . (19) 

As a result, the pseudospin operators take the form 

^1 = I (oIlCjL - cl^c.n] ■ (20) 



'j ~ 2 V ^'^^ ^'^^ 

The latter representation makes it clear the meaning of the pseudospin operators, showing 
that Sj describes the atomic tunneling between the wells of a double-well potential in the 
j-th lattice cite, S^ corresponds to the Josephson current between the wells, and the operator 
Sj defines the atomic imbalance between the left and right wells in the j'-the lattice cite. 
Employing these pseudospin operators yields [15] the single- atom energy part 

H, = Eo- QS^ , (21) 

in which Eq = {Ei +E2)/2 and Q !^ E2 — Ei is the interwell tunneling frequency. The binary 
interactions enter the part 

-"«i ~ 2 -^ ~'~ *■? « J ~ ^ij^i ^j ; (22) 

with Aij being the direct interaction between atoms in the i-th and j-th cites, Bij, the inter- 
actions of atoms in the process of their tunneling between the wells, and lij, the interactions 
of atoms shifted in one of the wells. All these quantities are explicitly expressed through the 
corresponding matrix elements of atomic interaction potentials over Wannier functions, as 
is done in Ref. [15]. 

For creating different states, it is important to have the possibility of varying atomic 
interactions in a wide range. In this regard, trapped atoms in optical lattices provide a 
unique opportunity. First of all, nowadays one is able to trap degenerate atomic clouds of a 
variety of atomic species possessing different scattering lengths. These are the whole series 
of bosonic atoms, ^H, ^He, ^Li, ^SNa, ^^K, ^^K, 52Cr, ^SRb, ^^Rb, ^^scs, i^oyb, i74Yb (see 
the summary in Refs. [1-5,38-44] and references therein) and also ^°Ca [45], ®*^Sr [46,47], 
and ^^Sr [48]. There are experiments with several degenerate trapped Fermi gases, '^He, 
^Li, *^°K, and ^^^Yb (see review articles [49-51]). It is also possible to create different atomic 
molecules, made either of bosons, ^^Na2, *^Rb2, ®^Rb2, ^^^Cs2, or fermions, ^Li2, ^°K2, as well 
as heteronuclear molecules, such as Bose-Fermi molecules *^Rb-'^°K or Bose-Bose molecules 
^^Rb-^^Rb and *^^K-*^Rb [52]. The strength of atomic interactions can be varied in a very wide 
range by means of the Feshbach resonance techniques [6,7,38,52,53]. Finally, there are atoms 
with long-range dipolar interactions, such as ^^Cr possessing large magnetic moments [54,55], 
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and also Rydberg atoms [56] and polar molecules [57]. Moreover, in optical lattices, it is 
admissible to regulate the shape of the double-well potential, varying by this the tunneling 
frequency [15,58]. Therefore, the system parameters can really be varied in the range of 
several orders of magnitude. 

Employing the pseudospin representation, discussed above, for the effective Hamiltonian 
(14), we obtain 

H, = w.NEo + yY1 ^^^- - ^-^ Y.^^ + "^'Yl {B^,StS'; - h.StS^) . (23) 

To some extent, this Hamiltonian is analogous to the effective Hamiltonian describing solids 
with pores and cracks [59] , in which the latter are the mesoscopic regions of disorder inside 
a regular solid matrix. 

For the insulating double-well lattice, the role of the order operator is played by the 
pseudospin operator 5*1, whose average gives the mean atomic imbalance. According to Sec. 
I, there are two kinds of averages, one defining the ordered state and another, disordered 
state. For an operator A, these two types of the averages are 

TrK.exp(-/3if^) 
where v = 1,2. For the order operator, one has 

(^l)i^O, (5|)2 = 0. (25) 

This means that the ordered state is characterized by a nonzero order parameter, while the 
disordered state, by the zero order parameter. 

4 Mixture Order Parameters 

It is convenient to define the reduced transverse averages 

a;, = 2(5f), , |/, = 2(5f),, (26) 

and the normalized order parameter 

su ^ 2{S^)^ . (27) 

To study the behavior of these averages under varying system parameters, let us resort 
to the mean-field approximation. Then, by using the standard mean-field techniques (see, 
e.g., [60]), for the free energy (12), we find 

F = F, + F2, F, = E,-NT\n r2cosh ^ , (28) 

where the first term in F^ is 

E, = w.EoN + ^[A-Bxl + Isl) N (29) 



and in the second term, the notations 

K = w, ^^l + wlPsl , Q, = Q- w,Bx, (30) 

are used. The parameters A, B, and / are given by the averages 

^-^E^^- ^-l^E^^- ^-^E^^- (31) 

*7^i *7^i *7^i 

For the mean tunnehng intensity, defined in Eq. (26), we get 

The average Josephson current is zero, 

yu = o, (33) 

as it should be in equihbrium. And for the order parameter (27), we find 

Su = wls^— tanhi^j . (34) 

Let us introduce the dimensionless tunnehng frequency and the tunnehng interaction 

parameter, respectively, 

^ , B 

Also, let us define the dimensionless quantities 

(^u = y-j^ = uj- bw^x^ , K = j-^ = Wy^bjl + (1 - hYwlsl , (36) 

corresponding to Eqs. (30). In what follows, let us measure temperature in units of / + 5. 
Then the tunneling intensity (32) reads as 

x. = ^. ^ tanh (^ (37) 



and the order parameter (34) becomes 



^ t-h (I 1 , (38) 



In view of the definition of the ordered and disordered states, symbolized in conditions 
(25), we have 

si 7^ , S2 = . (39) 

Therefore for the ordered state, from Eqs. (38) and (39), we get the tunneling intensity 

xi = — , (40) 



while the mean atomic imbalance si is defined by the equation 



2/. ,^. . I hi 



hi = wi{l-h)t^nh\^\ . (41) 

For the disordered state, for which the mean imbalance is zero, S2 = 0, the tunneling 
intensity is 

X, = tanh (^) . (42) 

In addition to the order parameters si and S2, the mixed system is characterized by the 
quantities wi and W2, playing the role of the geometric weights of the corresponding states. 
These weights, according to the theorem of Sec. 2, are defined as the minimizers of the free 
energy (12), under the normalization conditions (15). This minimization yields 



2u + wi^i — UJ2X2 2m — wiXi + ^22^2 ~ (1 ^ b)sl 

4u-{l-b)sl ' ^' " 4m - (1 - b)sl 



ZjU t u/ixi — 1^2-^2 ^" ~ uJiJbi T uJ2-i-2 ~ \^ ^ "j-'i / ,o\ 

^1 = —rr. — n — iJ72~ ' ^2 = — — r, — rT72 ' (^3) 



where the notation 

«~ (44) 

is used. The inequality (PF/dwf > results in the condition 

4:U + b{xj + xl) -{l-b)sj>0 . (45) 

The thermodynamics of the mixed system is described by the set of Eqs. (39) to (43). 

5 Quantum Phase Transition 

Setting T = in the above equations gives the tunneling intensities 

xi = — , X2 = l (46) 

Wi 

and the order parameters 



si = Wl- 4' ^2 = 0. (47) 

V ^1 

The state weights (43) become 

2u + b-uj 2u + b + uj -I ,^. 

"^ = 4. + 26-1 ' ""' = Au + 2b-l ■ ('') 

At zero temperature, the free energy coincides with the internal energy: 

F = {H) (T = 0) . (49) 

It is convenient to introduce the dimensionless internal energy 

^''"'' - WTW, • (5°' 
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for which we find 

1 ll]-\ 111 

E{wi) = eQ + -{2u + h-2uj-uj'^) -{2u ^h-uj) + -^{Au + 26 - 1) , (51) 

where 

The minimization procedure now imphes that 

^.0. ^!f^>0. (53) 

owi owl 

The first of these equations leads to Eqs. (48), while the second gives the inequality 

4m + 26 - 1 > . (54) 

By its definition as a matrix element over Wannier functions [15], the interaction strength 
B is assumed to be positive and less than I . Then, in view of definition (35), one has 

0<6<i (0<5</). (55) 

Taking into account condition (54) yields 

1 / l-26\ 

< Mo < - [u,^ ^^ j . (56) 

In order that Si in Eq. (47) be real, it is necessary that u be less than Wi. With wi from 
Eq. (48), we have 

UJ<\ (si > 0) . (57) 

Inequalities (54) to (57) specify the region of the system parameters, when mesoscopic dis- 
order can arise in a double-well optical lattice. 

For the double- well lattice with mesoscopic disorder, under the state weights (48), the 
internal energy (51) takes the form 

The energy of the mixed state should be compared with that of the pure states. If the system 
would be completely ordered, so that wi = 1 while W2 = 0, then the energy would be 

^(1) = eQ + -{2u + h-uj'^ -I) . (59) 

And, if the system would be absolutely disordered, so that Wi = while W2 = 1-, then we 
would have 

E(O) = eo + -[2u + h-2uj-uj'^) . (60) 
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From Eqs. (59) and (60), we get 



E(0) - i?(l) = ^-^ , (61) 



which shows that, under condition (57), the energy of the pure ordered state is lower than 
that of the pure disordered state, that is, the ordered state is preferable, being more stable. 
To compare the pure states with the mixed state, we consider the differences between 
the related energies, 

E{1) - E{wi) = {u - uo){l - wi)\ E{Q) - E{wi) = {u - uo)wl . (62) 

These equations tell us that the quantum phase transition could occur either at m = mq oi 
when wi = 1. From Eqs. (48) and (56), we see that the weight wi can be written as 

». = i + si^. (63) 

2 8(-u — Mo) 

Since, by definition (15), wi > and wi < 1, we find that 

U > Uq {wi > 0) , 

u>Uc {wi < 1) , (64) 

where 

Uc = uo^ — . (65) 

Under condition (57), Uc > Uq- Hence, the mixed state exists only if m > Uc- 

In this way, the quantum phase transition occurs at -u = Uc, given by Eq. (65), when 

Wi = 1. For u > Uc, the mixed state is more stable, while for u < Uc, the pure ordered state 

becomes preferable. When u — )■ Uc, then wi — )■ 1. And if -u — )■ oo, then wi — ?■ 1/2. 

The parameter u, defined in Eq. (44), has the meaning of the disorder parameter. For 

sufficiently large u > Uc, the regions of mesoscopic disorder appear in the system, making 

this mixed state preferable as compared to the pure ordered state. 

6 Temperature Phase Transition 

In the pure system without mesoscopic disorder, when wi = 1, the order parameter si, 
describing the mean imbalance between the wells of a double-well potential, tends to zero at 
the temperature 

Tc = ^^~^t , ■ (66) 

2arctanh(w) 

But, if the system is mixed, such that wi is described by Eq. (48), then the transition 
temperature, at which si — )■ 0, is 

^ (1 - h)uj 
' 4arctanh(2a;) ' ^ ' 

which is lower than temperature (66). Moreover, the situation for the mixed lattice is not 
as simple. The phase transition, at rising temperature, can happen not as a second-order 
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transition at T* but as a first-order transition at some intermediate temperature Tq between 
temperatures (66) and (67), provided that at this temperature Tq the mixed state becomes 
less favorable as compared to the disordered state. 

To study the behavior of the mixed state at finite temperatures, we take into account 
that the tunneling parameters u and h are small, in agreement with inequalities (55) and 
(57). Hence, for numerical calculations, these parameters can be neglected. Then Eqs. (66) 
and (67) reduce to 

Tc = \. T; = 1 . (68) 

For convenience, we define the weights as 

Wi = w , W2 = 1 — w . (69) 

And we solve numerically the equations for the order parameters 

s, = tanh ^-^j , w = ^^^-^ . (70) 

These equations, generally, can possess several branches of solutions. Among them, we have 
to choose that one corresponding to the smaller free energy. Also, we need to compare the 
free energy of the mixed state with that of the pure system, for which w = 1 and the order 
parameter is given by the equation 

s, = tanh (^) iw = l). (71) 

In addition, we have to consider the disordered state, for which si = and w — )■ 1/2. 
Comparing the free energies of all possible solutions, we select that solution whose free 
energy is minimal, which corresponds to the most stable system. 

In Figs. 3 and 4, we show the behavior of the stable branches of solutions as functions 
of temperature for different disorder parameters u. Figure 5 presents the dependence on the 
disorder parameter u of the order-disorder transition temperature Tq. The overall picture is 
as follows. 

(i) Nonpositive disorder parameter: 

M < . (72) 

Mesoscopic disorder does not appear at all. For temperatures below the critical T^, the pure 
ordered state is stable, 

w = l (0 < T < TJ . (73) 

At Tc, the second-order phase transition takes place from the ordered state, where si ^ 0, 
to the disordered state, where s\ = 0. 

(ii) In the region of the disorder parameter 

<U<Un , (74) 

where «„ ~ 0.45, the pure ordered state remains till the temperature Tq, 

w = l (0 < T < To) , (75) 
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where there happens a first-order phase transition to the disordered state. The transition 
temperature is in the interval 

T; < To < T, , (76) 

depending on the value of u, as is shown in Fig. 5. 

(iii) In the region of the parameter 



Un<U < - , (77) 



there exists the nudeation temperature 



Tn = ^^^ , (78) 

2arctanhv2-u 



below which the pure ordered state is stable, 

w = l (0 < T < T„) . (79) 

But, at Tn, there appears mesoscopic disorder, and the system becomes mixed, 

w<l {Tn<T<To). (80) 

At To, such that 

T„ < To < T, , (81) 

there occurs a first-order phase transition from the mixed state to the disordered state. 
When u — > «„, then T„ — )■ Tq. 

(iv) For the parameter 

u = l, (82) 

the nucleation temperature T„ = 0. At zero temperature, the system is purely ordered, 

w = l (T = 0) . (83) 

But it is mixed for all temperatures above zero up to To, 

w <l (0 < T < To) . (84) 

At To, there is a first-order phase transition from the mixed state to the disordered state, 
(v) In the region 

-<u<- , (85) 

the system is always mixed, 

w <1 (0 < T < To) , (86) 

up to the temperature Tq, where a first-order phase transition happens. 

(vi) For the value 

u = l, (87) 
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the system is mixed up to the temperature T*, 

w <l (0 < T < T;) . (88) 

The temperature T* corresponds to a tricritical point that is the end-point on the hne of 
first-order phase transitions, separating these from the hne of second-order transitions. 

(vii) At large values of the disorder parameter, such that 

u>\. (89) 

the system is always mixed, 

w <l (0 < T < T;) . (90) 

At temperature T^*, there occurs a second-order phase transition to the disordered state, 
where si = 0. 

This analysis shows that the regions of mesoscopic disorder arise in a double-well lattice 
when the disorder parameter u is sufficiently large. By its definition, the parameter u is 
the ratio of the disordering interactions to the effectively ordering ones. The appearance 
of the mesoscopic disorder is due to the competition between the ordering and disordering 
interactions. 

7 Conclusion 

The double-well optical lattice in the insulating state has been considered. In such a lattice, 
depending on the system parameters, there can exist the ordered state, when the mean 
atomic imbalance is nonzero and the disordered state, with the zero atomic imbalance. 

We have shown, that, in addition to these pure states, there can arise an intermediate, 
mixed, state, when, inside a generally ordered lattice, there develop the regions of mesoscopic 
disorder. These regions are mesoscopic in the sense of their typical sizes being intermediate 
between the lattice spacing and the total system size. The shapes and locations of the 
mesoscopic regions of disorder are random, which requires to invoke the averaging over all 
their random spatial configurations. Accomplishing this averaging allows one to define an 
effective renormalized Hamiltonian and to calculate all thermodynamic characteristics. All 
mathematical details, describing the configuration averaging have been given in the earlier 
papers [21-28,61,62]. 

The considered mesoscopic disorder appears in a self-organized way and may happen 
without any external random forces. The stability of this kind of disorder is due to the 
existence in the system of different types of pair interactions, among which there should exist 
the interactions that do not favor the ordering. The ratio of the competing disordering to 
ordering interactions plays the role of a disorder parameter. The latter has to be sufficiently 
large in order to stabilize the mixed state with the mesoscopic disorder. 

As is discussed in Sec. 3, the optical lattices of trapped atoms provide a unique opportu- 
nity of regulating practically all system parameters. This can be done by varying the types 
of atoms, by creating atomic molecules, and by changing the lattice parameters. It is also 
possible to regulate the shape of the double-well potential, varying by this the tunneling 
frequency [15,58]. The strength of atomic interactions can be varied in a very wide range by 
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means of the Feshbach resonance techniques [6,7,38,52,53]. Some atoms possess interactions 
with long-range dipolar forces, such as ^^Cr [54,55] and Rydberg atoms [56]. There are also 
polar molecules with long-range forces [57]. It is admissible to create molecules of large 
size, whose interactions also have long-range nature, such as the interactions in liquids and 
plasmas [63], having, for instance, the Lennard- Jones form [64]. 

One more possibility of influencing the system states is by invoking external fields. Pump- 
ing energy into the system can produce effects that are similar to increasing the system 
temperature. It is known [65,66] that sometimes even weak external temporal perturbations 
essentially change the system behavior. The stationary action of alternating external fields 
can be described as the effective rise of temperature, which can lead to the appearance of 
mixed states, as in the statistical approach to quantum turbulence [67]. It may also be 
possible to excite a part of atoms out of the bound states and to create different asymmetric 
localized states above the double- well lattice, such as exist in bichromatice lattices [68]. The 
localization should occur in the disordered spacial regions. 

In order to give an understanding of how the effective temperature can be estimated, in 
the case of a stationary action of an external alternating potential, let us give an example. 
Suppose the system is subject to the action of an external modulating potential V(t), with 
amplitude Vmod, alternating with a period tmod during the time At. If the modulation period 
is much longer than the local equilibration time, then the system is in quasi-equilibrium at 
each instant of time. If there exists a thermostate with temperature T and, in addition, the 
modulating potential pumps energy into the system, then the effective temperature can be 
estimated as 

-^ eff -L -\- 1 fYiod ; 

where the modulation temperature is 

rri ^'' Tr ^mod a , 

-^ mod 7 ''rnod ~~Z ^mod^'^ ■ 

'^mod ^'^ 

The observation of the arising mesoscopic disorder and the level of disordering can be 
realized by optical measurements [69] or scattering experiments [70]. 

In this way, there are many possibilities to vary the system parameters. Therefore, it 
seems feasible to find the conditions favoring the appearance of the mesoscopic disorder. 
The level of disorder can also be regulated, which can be used in applications. 
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Figure Captions 



Fig. 1 Scheme of the ordered state in an insulating double-well optical lattice. 



Fig. 2 Scheme of the intermediate mixed state, with mesoscopic disordered regions inside 
the ordered matrix denoted by arrows. 



Fig. 3 Imbalance order parameter as a function of temperature (in dimensionless units), 
with the marked points of the order- disorder phase transitions, for different disorder param- 
eters: (1) u = 0; {2)u = 0.1; (3) u = 0.3; (4) u = 0.4; (5) u = 0.51; (6) u = 0.75; (7) u = I; 
(8) u= 1.5. 



Fig. 4 Geometric weight of the ordered phase in the mixed system as a function of 
temperature (in dimensionless units), with the marked points of the order-disorder phase 
transitions, for the same disorder parameters as in Fig. 3. 



Fig. 5 Order-disorder transition temperature (in dimensionless units) as a function of 
the disorder parameter u. The phase transition is of second order for -u = and u > 3/2; it 
is of first order for < m < 3/2; and u = 3/2 corresponds to a tricritical point. 
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Figure 1: Scheme of the ordered state in an insulating double- well optical lattice. 
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Figure 2: Scheme of the intermediate mixed state, with mesoscopic disordered regions inside 
the ordered matrix denoted by arrows. 
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Figure 3: Imbalance order parameter as a function of temperature (in dimensionless units), 
with the marked points of the order- disorder phase transitions, for different disorder param- 
eters: (1) u = 0; {2)u = 0.1; (3) u = 0.3; (4) u = 0.4; (5) u = 0.51; (6) u = 0.75; (7) u = 1; 
(8) u = 1.5. 
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Figure 4: Geometric weight of the ordered phase in the mixed system as a function of 
temperature (in dimensionless units), with the marked points of the order-disorder phase 
transitions, for the same disorder parameters as in Fig. 3. 
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Figure 5: Order-disorder transition temperature (in dimensionless units) as a function of the 
disorder parameter u. The phase transition is of second order for u = and u > 3/2; it is of 
first order for < m < 3/2; and u = 3/2 corresponds to a tricritical point. 
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